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O ■ Abstract 

o 



Multibubble solutions for a cosmological model which lead to ther- 
mal inflationary states due to a semi-classical tunneling of gravity are 
calculated. 



On 

"o 

W)' I. One domain solutions 

In earlier work done by several different approaches , an infla- 
tionary solution was found as semiclassical tunneling using a simple model 
of quantum matter and semiclassical gravity. The model comprised massless 
noninteracting scalar bosons conformally coupled to gravity in a RW universe 
with k = +1 and a positive phenomenological cosmological constant A. For 
such model there is a classically forbidden region connecting two classically 
allowed regions. The inner one of these two is inappropriate for our model 
since the inner region necessarily involves quantum gravity. The model was 
evaluated by two different approaches, first by an explicit statistical mechan- 
ics approach and then by a wave function of the universe approach. In both 
of these approaches the initial method was to calculate the solutions in the 
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outer classically allowed region. The result of tunneling is to produce an 
inflationary thermal solution in the outer region. This tunneling yields both 
a definition of physical time and entropy appropriate to an inflationary solu- 
tion. 

We use a RW metric written in conformal time coordinates: 

ds 2 = a 2 (u)[du 2 - dx 2 - sin 2 x<^ 2 ] = a 2 (u)^ ^ u dx^ dx v (1) 

The action obtained after going over to a Euclideanized metric as 
appropriate for a tunneling region is 

1 = I dxA[ h {h2 +a2 ~ A/3fl4) + - 02) ] = Ig+ Im (2) 

after carrying out a conformal transformation g = a 2r y and = aip, whence 
R/6 = 1. One then evaluated the norm of the wave function of the universe 
and showed that 

P(a) = [ dSVU, a) | 2 = [ D[a] [ £U] e -'(<^>^o) (3) 

J J 2ni Jwkb J 

where and a are periodic functions of f3' . The f3' integral fixes the total 
energy to be zero and its saddle point value is the reciprocal temperature 
when we have exchanged the order of integration as in the second part of 
equation (3). The saddle point evaluation of a involves replacing a in the 
integrand by its saddle point value a. The functional integral gives 

J £>[0]e-W) = e-^W (4) 

where F((3) is the Helmholtz free energy. Then P(a) = e s = e S(3+SM , where 
S is the total entropy and Sq and Sm are respectively the entropies of the 
gravity and of the matter. The common saddle point solutions of the a and 
(3' integrations leads to the equation 

p = 2 r dx i ^ 2 4 

where x = (| ; ) 1 / 2 a and e is the thermal average energy multiplied by ^; 

e = §7rAA; B /r 4 = af3~ A . 
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Figure 1: Solutions of equation (5) 



Seeking the self consistent solution in the tunneling region one finds two 
solutions in a certain range of the values of x which can most clearly be 
seen by the following graph, where the dashed curve is the locus of solutions. 
Each point of the dashed line is the end point of a straight line from x^ 
at that energy. Thus we see of all there are no solutions for x less than 
some minimum value Xq and that there are two branches of the solutions for 
Xq < x < Xm = 1/ v^2- Both branches correspond to an increasing entropy as 
a function of x and both are maximum entropy states of nearly equal weight 
near the maximum value of the upper solution. The upper solution joins 
the real time domain where it corresponds to a static Einstein solution. The 
lower one emerges at higher entropy and corresponds to a thermal inflation- 
ary solution. Note that in this case there is no need for time to establish the 
thermal solutions, since it begins in thermal equilibrium when it emerges in 
the real time domain. 

In previous work we discarded the upper solution since it is of lower 
entropy and is also unstable. Subsequently we began to consider whether it 
is justified to ignore alternative saddle point solutions which cannot imme- 
diately be discarded as being of much less weight. The present paper is an 
effort to consider another alternative. 
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II. Multidomain solutions 



In the standard functional integral approach, if there is more than 
one saddle point the wave function in lowest order is a superposition of the 
several contributions. This method is not viable in our case since we are 
considering thermally mixed states. We apply instead another approach in 
which the superposition consists of the different solutions cooexisting in spa- 
tially separated regions. We worked out in detail the case of two domains; 
the multibubble solution were treated only approximately. Let us first con- 
sider the two bubble case. The two states are a Robertson- Walker (RW) 
inflationary universe and a static Einstein universe (SE). The two different 
solutions are matched on a singular surface, which is taken to be spherical. 
The Euclidean line elements in the two regions take the form: 

ds RW = \ x2 ( M ) i du2 + d XRW + sin2 XRwdQ 2 ] (6) 

ds 2 SE = ^[N 2 SE {u)du 2 + d X 2 SE + sin 2 Xs E dQ 2 ] (7) 

The singular hypersurface is given either by the equation xrw — 
fnw(u), or by the equation Xse = fsE^u)- The matching conditions on the 
hypersurface lead to the following relations: 

x(u) sin f RW (u) = -y= sin f SE (u) (8) 

x 2 (u)(l + f 2 RW (u)) = \{N 2 SE {u) + f 2 SE {u)) (9) 

Another constraint follows from the junction conditions that relate 
the jump in the Einstein tensor to the surface energy 0. As a toy model 
we assume that the surface energy is purely dynamic and it comes from the 
surface term in the action. We showed elsewhere || that the variation of the 
action with respect to the induced metric on the hypersurface yields, in our 
case the following constraint: 

fnw — constant. (10) 

The thermodynamic entropy we obtained is the contribution from the 

two states 

S = Srw + Sse (11) 
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each contribution being the sum of the gravitational and the matter entropies. 

rx+(e) _ ^_ _ 4 

lx-{e) 

1 rP , s 1 
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Jx-(e) 3 



(12) 



S SE (e) = 2A 



, rfwiV 2 (a;)^(x) + ^\/ 2 (x + )a/3 2 - 3 (e,/5) 
4 jo 3 



(13) 



Here Vrw = Jo" w sin 2 x dx and Vg^ = ff SE r u ) srn2 X are the volumes 
occupied by the RW and SE bubbles respectively. The function x(u) is the 
saddle point solution of the variational equations for the RW region: 

x 2 = x 2 - x 4 - e (14) 

with e being found from the equation 

r°)*= 2 r w , dx ' as) 

\ej J x _(e) ^/x' 2 -x'^-e y ' 

The (proper) temperature in the SE region is related to j3 with the lapse 
function: 

(3 2 (e,[3) = [ duN 2 {x). (16) 
Jo 

The condition for the two bubbles to be in thermal equilibrium is 

&(e,(a/e)*) = (4a)* (17) 

Equations (8)-(10) , (15)— (17) constitute the conditions to be satis- 
fied in equilibrium. Having found e, fnw, x(u), Jse{u) and N S e(u) from these 
equations, we can find the value of the entropy at equilibrium and check the 
stability from eqs. (11)— (13). Analysis of eqs. (15)— (17) shows that an equi- 
librium can only exist for the range of parameter a, it 4 /A < a < 7r 4 . In this 
range there is a critical value a c , such that the equilibrium is stable against 
global fluctuations when a < a c and is quasistable when a > a c . In either 
case the entropy of the two bubble state is greater than the entropy of the 
single RW universe. 

We now turn to the case of multiple bubbles. There are two possibili- 
ties: many RW bubbles in a SE matrix or many SE bubbles in a RW matrix. 
It can be shown that conditions for equilibrium remain the same as in the 
case of two bubbles. In our simple model the number of bubbles cannot ex- 
ceed 25. Moreover, the constraints limit the size of the RW region, so that 
inflation cannot last long. We are currently exploring possible extensions of 
the model to overcome these problems. 
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